Electronic Transport in Metal Nanocrystal Arrays: The Effect of Structural Disorder 

on Scaling Behavior 
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We investigate tlie impact of structural disorder on electronic transport in gold nanocrystal mono- 
layers. Arrays ranging from void-filled networks to well-ordered superlattices show clear voltage 
thresholds (Vt) due to Coulomb blockade, and temperature-independent conduction indicative of 
quantum tunneling. Current-voltage characteristics of arrays with and without long-range struc- 
tural order were found to collapse onto distinct scaling curves. The former follow a single power-law: 
I {V — Vt)^, C = 2.25 ± 0.1. The latter show additional structure, reflecting the underlying dis- 
ordered topology. 



The intriguing electronic and optical properties of in- 
dividual nanocrystal quantum dots have unleashed a 
flood of interest However, despite such discoveries 
as metal-insulator-type transitions in squeezed nanocrys- 
tal monolayers |^ and spin-dependent transport in mag- 
netic particle assemblies the simplest nanoparticle 
array — a single layer of metal nanocrystals — has re- 
mained poorly understood. The main reason is that the 
transport characteristics are strongly affected by three 
types of disorder: global structural disorder in the ar- 
ray topology, local structural disorder in the interparti- 
cle couplings, and local charge disorder due to random, 
immobile charges in the underlying substrate. Theoret- 
ical approaches investigating tunneling transport so far 
have concentrated on local charge disorder only A 
full treatment of the combined types of disorder is not 
available, even though large differences between spatially 
ordered and disordered structures might be expected due 
to the sensitivity of percolative charge transport phenom- 
ena to array topology and the exponential dependence 
of local tunneling resistances on the interparticle spac- 
ings. With nanocrystal arrays as "artificial solids" ||] 
expected to provide useful analogues and tunable test- 
beds for various bulk correlated electron systems, a full 
understanding of metal nanocrystal monolayers is vital. 

Experimentally, a high degree of structural order has 
been elusive for arrays between in-plane electrodes, and 
previous investigations of electrical conduction in 2D 
nanocrystal systems have been performed only on small, 
highly disordered, or multi-layered arrays Using 
newly-developed self-assembly techniques, we have fabri- 
cated large highly-ordered monolayers of dodecanethiol- 
ligated gold nanocrystals on substrates with in-plane 
electrodes. Transport measurements and subsequent 
transmission electron microscopy (TEM) on the same ar- 
rays allowed for direct correlation of electronic and struc- 
tural characteristics. Comparison of these superlattices 
with void-filled networks for the first time delineates the 
roles played by the different types of disorder. 

Nanocrystals were deposited on silicon substrates 
coated with lOOnm amorphous silicon nitride (Si3N4) 
(Fig. la). Under a 70/xm x 70/im area, the Si was 



etched away to leave a freestanding Si3N4 membrane 
"window," allowing TEM imaging Q. Thin (20nm) Cr 
electrodes were patterned using electron-beam lithogra- 
phy. 1-dodecanethiol ligated gold nanocrystals were syn- 
thesized as described in Ref. and dissolved in toluene 
to a concentration of about 2.4 x 10^'^ ml~^. The gold 
core radii varied from sample to sample (2.2-2.9nm), but 
for each sample were monodisperse to within 5%. 

Two different array preparation techniques yielded two 
distinct classes of arrays: samples with and without 
large-scale structural disorder. Simple deposition of 15- 
20 fil of colloid onto a substrate produced, upon drying, a 
nanocrystal monolayer of well-packed regions with short- 
range order, coexisting with numerous voids (area frac- 
tion 15-20%) throughout the 2D plane (Fig. lb). Higher 
particle concentrations produced an increased amount of 
multi-layered regions without increasing the long-range 
order. However, addition of excess dodecanethiol (vol- 
ume fraction 6.3 x 10~^) to the solution before deposition 
increased the nanocrystal mobility on the Si3N4 surface 
and prevented rapid dewetting of the solvent from the 
substrate. As a result, arrays with significantly smaller 
void or double layer fraction (about 5% combined) and 
excellent long-range periodicity could be self-assembled 
(Figs. lc,d) [|ll[. Electronic properties of a total of 14 
arrays were measured (7 prepared without and 7 with 
excess ligand), with dimensions defined by the electrode 
separations (200nm < D < 700nm) and widths (fixed at 
2fim). The resulting N x M arrays ranged in length from 
iV = 30 to 90 particles and were M w 270 nanocrystals 
wide. 

The samples were cooled, in vacuum, to below 77K to 
avoid parasitic conduction through the substrates. DC 
current-voltage (IV) characteristics were measured using 
Keithley 614 electrometers and a voltage source. Control 
measurements on substrates without monolayers showed 
no detectable currents (< 0.02pA) up to ±20V. We im- 
aged each array by TEM after the transport measure- 
ments. From analysis of the center-to-center distances, d, 
and the particle radii, r, distributions of the interparticle 
spacings 2s = d — 2r (Fig. la) were obtained, resulting in 
s = 0.85 ± O.lnm for the arrays with voids. For the well- 
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FIG. 1. (a) Sketch of a nanocrystal monolayer and 
in-plane electrodes (not to scale), and of the interparticle ge- 
ometry, (b) TEM image (detail) of a typical array formed 
Without excess dodecanethiol, showing many voids and an ab- 
sence of long-range order. The 2D Fourier Transform (2DFT) 
is of a larger area, (c) Typical superlattice formed with ex- 
cess dodecanethiol, showing <5% voids and long-range or- 
dering, (d) Highly ordered superlattice (and 2DFT) between 
electrodes visible at the upper left and lower right. The ran- 
dom nanocrystal image intensities are due to random Bragg 
diffraction of the electron beam. 

ordered superlattices the excess dodecanethiol increased 
the mean spacing to s = 1.2 ± O.lnm. 

The IV curves (Fig. 2) were highly symmetric, of the 
same overall shape for all samples, and showed no hys- 
teresis at the slow ramp rates used (<5 mV/s). The 
strongly non-ohmic behavior, in which current flow re- 
quires the applied voltage to exceed a finite thresh- 
old, is characteristic of Coulomb blockade of transport 
P,[7|JT^-|r^ . The electrostatic energy needed to add one 
electron to a quantum dot of charge q and (self-) capac- 
itance Co = ATTeeor leads to a single-particle Coulomb 
blockade voltage, Vq = q/Co, below which tunneling is 
suppressed (Fig. 2, lower inset). For V > Vq, cur- 
rent flows with resistance R — dV/dl; R ^ h/e^ for 
our nanocrystals. The overall IV characteristic then 
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FIG. 2. IV curves for a typical superlattice {D = 330nm 
X width 2/im). The upper inset magnifies the data. The lower 
inset shows a schematic IV curve for a single nanocrystal, as 
described in the text. 

arises from the series-parallel combination of many tun- 
neling paths from particle to particle throughout the ar- 
ray. Random, parasitic charges in the substrate — local 
charge disorder — induce the effective charges q S (0, e) 
on the nanocrystals, randomly placing the local Vq in 
the interval (0, e/Co), leaving R unchanged. Local struc- 
tural disorder, on the other hand, produces variations 
in the tunnel distance, 2s, and thus a wide (most likely 
log-normal) distribution of i?, leaving Vb essentially un- 
changed ||l5|. The exponential dependence of i? on s 
makes a large variance in the interparticle resistances al- 
most a certainty, even in well-assembled arrays (e.g. Fig. 
Id). 

Using the dielectric constant e « 2 for dodecanethiol, 
and q = e/2, we find Co ~ 0.5aF and Vq « 150 
mV for a typical Au nanocrystal. The capacitance be- 
tween neighboring particles iQ, based on the geome- 
try observed by TEM, is C12 = 0.25aF < Cq. In this 
regime of small interparticle capacitive coupling, hys- 
teresis due to long-range charge re-arrangements is not 
expected [Q. Furthermore, for all accessible tempera- 
tures e^/ max{Co, C12} » keT and thus the array trans- 
port properties are essentially temperature-independent 
(Fig. 2). This behavior implicates direct, interparticle 
quantum tunneling as the conduction mechanism and is 
in contrast to the strong temperature dependence ob- 
served in other nanocrystal systems with large Coulomb 
blockade voltages [^jj^. Arrhenius behavior in those sys- 
tems we believe is due to activated hopping via trapping 
sites provided by electron-rich 7r-bonded linker molecules 
attached to the metal particle cores. This difference 
between tunneling and activated hopping suggests that 
chemical modifications to nano crystalline systems may 
change fundamentally the physics governing transport, 
and highlights two different routes towards control of 
nanostructure conductance: physical tailoring of sepa- 
rations between metal cores and chemical manipulation 
of molecular links between the cores. 

How charge disorder alone affects the global IV charac- 
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teristics measured across large arrays has been calculated 
by Middleton and Wingreen (MW) 0. They find a well- 
defined, global threshold voltage, Vt, which delineates a 
second order phase transition: for V < Vt the array cur- 
rent is zero, and for V > Vt the array conducts with a 
power-law scaling form, / ^ M{Vo/R)[{V — Vt)/Vt]^ . In 
Figure 3 we plot the normalized current, i = IR/{MVo), 
as a function of normalized voltage above threshold, 
V ~ {V — Vt)/Vt, for all 14 measured arrays. The ar- 
ray width M was obtained from the TEM images. R 
we adjusted to obtain collapse of the curves. Typical 
R were around 50 and 300Tri for arrays without and 
with excess dodecanethiol, respectively, the increase cor- 
responding to the increase in average particle spacing de- 
scribed above. These R are in accord with conducting-tip 
atomic force microscope studies of tunneling through self- 
assembled alkanethiol monolayers |17|. Consistent choice 
of Vt (not far from the values estimated by eye from the 
IV curves, e.g. Fig. 2) produced a robust scaling behav- 
ior of the IVs for the ordered arrays assembled with ex- 
cess ligands (Fig. 3a), yielding an average C = 2.25±0.1. 

The conductivity exponent C is related to the mean- 
dering of current paths in the charge-disordered land- 
scape. Even in regular 2D arrays current paths are not 
straight, but exhibit transverse fluctuations that extend a 
distance ^± oc u^'' due to the quenched disorder |^. The 
total current above threshold across an array of width 
M is proportional to the number of independent parallel 
paths, M/^x, which leads to i oc vM/£^± oc Mv'^'^^, i.e. 

— T] + 1. The determination of the transverse correla- 
tion length exponent rj arises in many physical phenom- 
ena connected to interface growth or directed percolation. 
MW argue that rj = 1/z, where z = 3/2 is the roughness 
exponent for Kardar-Parisi-Zhang (KPZ) models for in- 
terface growth in 2D ; thus C ~ 5/3 . Simulations 
show C <« 2 (i.e. ri 1) for square arrays of up to 
400 X 400 particles Q. The power-law form of the data 
in Fig. 3a suggests that scaling models are indeed valid, 
but the measured C is significantly larger than predicted 
by available theory. 

Several previous experiments have reported a wide, 
sample-dependent spread in the scaling exponent: 1.4 < 
C < 2.0 for lithographically patterned junction arrays 
[^,|3|, 1.6 < C < 2.1 for polydisperse nanocrystals [||, 
and 2.2 < ( < 2.7 for small multi-layered arrays This 
contrasts with the highly reproducible values we observe 
for C in the structurally ordered arrays. However, we 
find that arrays with large void fraction display i{v) scal- 
ing behavior with characteristic slope changes as well as 
larger sample-to-sample variations (Fig. 3b). The differ- 
ences between the two classes of IVs may be understood 
in terms of changes in network topology. At large void 
fraction, neighboring voids produce bottlenecks, locally 
cutting off the transverse correlation length In the 
extreme case, conduction is reduced to several parallel 
ID channels, each with linear i oc v but threshold Vt dis- 



tributed over some range. For V barely larger than the 
smallest Vt only one channel is open and the overall IV 
is linear. As the applied voltage is increased, there will 
be a cross-over region in which a growing number of par- 
allel channels conduct. This is born out by the simulated 
IV characteristic of such a (l-l-l)D system (Fig. 3b, in- 
set). Depending on the number of independent channels, 
behavior resembling power-laws in the cross-over region 
with exponents ranging from CR;ltoC>2.5 can easily 
be reproduced [Q. For sufficiently large v all ID chan- 
nels have opened and the overall IV necessarily becomes 
linear again. Actual arrays most likely are amalgams of 
locally 2D patches connected by ID channels. Therefore, 
once all ID bottlenecks are filled, the overall IV charac- 
teristics are dominated by the remaining 2D patches and 
C approaches a value close to that of ordered 2D arrays. 
This is the behavior seen in Fig. 3b, where ( « 2.7 in 
the cross-over region turns over to ^ = 2.16 ± 0.1 beyond 
w « 2. 

Even for our largest accessible v, neither the IVs of the 
ordered nor disordered arrays turn linear (Fig. 3a, b), 
as would be expected once approaches a single lat- 
tice spacing. This finding is consistent with results from 
lithographically patterned arrays below u w 10 and 
shows a remarkable extent of the scaling regime. 

The thresholds Vt obtained from the scaling collapse 
of the ordered arrays grow linearly with array length, N 
(Fig. 3c): VT/{e/Ca) = aN with a = 0.25 ± 0.02. Sim- 
ulations |4| show that a for a given lattice depends only 
on the capacitive coupling between neighboring particles 
and decreases as coupling increases; 2D square arrays in 
the limit C12/C0 give a = 0.338. Due to both the 
larger coordination number in our hexagonal arrays and 
the finite C12/C0 « 0.4 we expect a < 0.338, consis- 
tent with the measured value. For arrays with large void 
fraction, N is poorly defined, leading to strongly sample- 
dependent Vt(A^) values. 

From our findings, two key results emerge. First, suf- 
ficiently large amounts of topological disorder, due to 
voids in the monolayer, lead to distinct deviations from 
simple power law behavior in the IV characteristics. In 
situations where direct imaging is impossible, detailed 
examination of /Vs, therefore, may provide clues about 
an array's large scale topology. Second, IV characteris- 
tics of long-range ordered arrays are well-fit by a single 
power law, despite the existence of both charge disorder 
and an inherent wide distribution of tunnel resistances. 
We believe this indicates that while R may be exponen- 
tially sensitive to variations in the interparticle separa- 
tion, charge disorder nevertheless plays the dominant role 
in selecting optimal current paths across the array. The 
reason most likely lies in the extremely non-linear, local 
IV characteristics (Fig. 2, inset), which effectively shut 
off all current flow unless Vq is exceeded . Within this 
picture, our finding of an exponent C w 2.25 in the struc- 
turally well-ordered arrays implies rj ~ 1.15 and thus, as 
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FIG. 3. Scaling behavior of IV curves. (a,b) Log-log plots of all data with normalized voltage beyond threshold 
V = {V — Vt)/Vt and current i = IR/{MVo) as described in the text, (a) Data from seven monolayers with long-range 
structural order. The solid line shows best-fit power-law i — v'' with ( = 2.25. (b) Data from seven disordered monolayers. 
The solid line shows = 2.25 as in (a). Inset: IV curve from a (H-l)D simulation of 100 independent parallel channels each 
100 particles in length, (c) Threshold voltage Vt, in units of e/Co, versus array length N. The line is a fit to Vr ~ aN{e/Co), 
which gives a = 0.25 for the superlattices formed with excess dodecanethiol. 



Vt is approached from above, a stronger divergence of 
than would be expected from charge disorder alone 
(where r] « 0.67). In other words, the spread in tunnel 
resistances appears to produce a more rapid growth of 
transverse fluctuations in the meandering current paths. 
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